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Abstract 

The structure of the gaussian auxihary field approximation in the 
theory of phase ordering kinetics is analysed with the aim of plac- 
ing the method within the context of a systematic theory. While 
we are unable to do this for systems with a scalar order parame- 
ter, where the approximation remains uncontrolled, a systematic 
development about the gaussian approximation can be outlined 
for systems with a vector order parameter in terms of a suitably 
defined 1/iV-expansion. 
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1 - Introduction 



Although much progress has been made in the understanding of phase order- 
ing kinetics^'^, from the point of view of the theorist basically this remains an 
unsolved problem. The reason is that in the case of a scalar order parameter, 
which is the most relevant for experiments, no systematic scheme for the de- 
velopment of a perturbation theory is available. This requires the existence of 
a soluble zero order approximation which accounts, at least qualitatively, for 
the relevant physical features of the problem and of a well defined procedure 
for the calculation, at least in principle, of the higher order corrections. The 
case of a system with continous symmetry is in better shape since the 1/N- 
expansion meets these requirements, at least in the case of a non conserved 
order parameter. For a conserved order parameter there are indications that 
the large-N limit might be singular^. 

Despite this very unsatisfactory situation, recently much progress has 
been made in the development of analytical methods for the computation of 
the structure factor, through extensive use of approaches based on the intro- 
duction of a gaussian auxiliary field"^'^ (GAF) which improve on the original 
idea of Ohta, Jasnow and Kawasaki^ (OJK). An exaustive critical account of 
these theories has been given by Yeung, Oono and Shinozaki^. The success 
of this approach, which for the moment is mostly limited to non conserved 
order parameter, amounts to the very accurate reproduction of the scaling 
function for scalar order parameter^'^'^ as known from experiments or nu- 
merical simulations, and to the prediction of power law tails in the case of 
vector order paramctcr^'^. Thus, these theories seem to incorporate those ba- 
sic ingredients that a real theory of phase ordering kinetics should have. The 
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shortcoming is that the assumption that the auxihary field obeys gaussian 
statistics is totally uncontrolled. More than theories for the moment these 
are sophisticated computational prescriptions which are justified a posteriori. 

A substantial progress toward a systematic theory then would be made 
if it were possible to identify a scheme within which a GAF approximation 
plays the role of the zero order approximation together with the expansion 
parameter which generates the higher order corrections. This type of project 
is illustrated in the paper of Bray and Humayun^. 

Motivated by these considerations here we overview the GAF approxima- 
tions with the aim of exposing those features which help to put into focus 
what is required for the eventual development of a systematic theory. For 
pedagogical reasons we begin with a detailed discussion of the one parti- 
cle problem which is exactly soluble and therefore allows to illustrate clearly 
what is involved in a GAF-type approximation. The same pattern of analysis 
than will be apphed to the field theory case. 



2 - One Particle 



Let us consider one particle in a double well potential and in contact with 
a thermal reservoir. The decay process from the instability point of this 
system has been thouroughly studied in the literature-^'^. In the limit of zero 
temperature the equation of motion for the position (pit) is given by 

= r0 - (2.1) 

with r > , g > 0. In order to study the quench of this system from high 
temperature to zero temperature, let us consider a gaussian probability dis- 
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tribution for the initial condition 0o = 4>{to) 



1 



Due to the symmetry of the problem the average position of the particle 
vanishes identically < >= and we concentrate on the behaviour of the 
fluctuations 

/ + OC 
d(t)P{(t>,t)(t)^ (2.3) 
-oo 

where P(0, t) is the probability that the particle occupies the position (j) at 
the time t. This quantity can be computed exactly since the equation of 
motion (2.1) can be solved 

with T = Thus, for the probabihty density we find 

df-\t-toA) 



p{(t>,t) = Po{r\t-to,(t>))- 



6Xp{— 2^ [l-(g/r)02{l-T-2)]} 



(2.5) 



v/2^fA^ [1 - {g/rW{l - T-2)]3/2 
where f~^{t — to, 0o) is the inverse of f{t — to, 0o)- Inserting into (2.3) 



r 



'^(^) = ^p3l) ~ V^e^il - erfiV^)]} (2.6) 

where x = r/[2Ag{T'^ — 1)] and erf{z) is the error function. For short time 
(2.6) yields exponential growth of the fluctuations 

S{t) ~ At' (2.7) 

due to the initial instability, while for large time we get 

S{t)^'-j^^^[l-V^ (2.8) 



which describes the saturation toward the finite equihbrium value S{oo) = 
^eq — ^/s* to the fact that eventually the particle sits at the bottom of 
one of the two potential wells and the probability density (2.5) develops two 



If the exact solution of the problem was not available, this type of be- 
haviour could not have been obtained via a straightforward perturbation 
expansion in the nonlinear coupling g. Zero order amounts to take the gaus- 
sian approximation in (2.5) which describes only the regime of exponential 
growth (2.7). Hence, the zero order theory does not reproduce the qualitative 
picture of the process, neither any improvement is obtained by taking into 
account corrections of finite order. The saturation to a finite final equilib- 
rium value is obtained within a perturbative scheme, as shown by Suzuki^°, 
by resorting to the infinite resummation of the most divergent terms in the 
series. 

However, rather than following this route, let us use the method which 
in the following will be generahzed to the field theory case. The idea is to 
introduce a new auxiliary variable m{t) through a transformation 



which takes care of the basic non-linear features of the problem in such a way 
that the behaviour of m(i) can be treated by straightforward perturbation 
theory. Namely, the transformation must be such that while m(i) is allowed 
to grow indefinitely the saturation of 0(t) to a finite value is induced by a. 
Substituting (2.9) into (2.1) one obtains the equation of motion for m{t) 



narrow peaks centered about the equilibrium values (peg — ±Jr/g. 




(f>{t) = a{m{t)) 



(2.9) 




(2.10) 
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with the transformed probabihty density of initial conditions 



Pmoimo) = Po{a{mo)\R)a'{mo) (2.11) 

where mo = m{to) and Po(0o|-R) is the probabihty density (2.2) conditioned 
to 00 belonging to the range of values R for which (2.9) is invertible. Thus, 
in terms of m{t) we cannot quite get the fluctuations (2.3), but fluctuations 
conditioned to (p E R 

S{t\R) = j d(t)P{(t),t\R)(t)'^ = j dmPm{m,t)a^{m) (2.12) 

where R is the domain 

R^ {<!>'< r/g) (2.13) 

and PmifTi^t) is the probability density of m at the time t. How important 
this restriction is depends on what is the statistical weight of trajectories 
lying outside R and this in turn is related to the size of the variance A of 
the initial probability density (2.2) compared to the size r/g of the domain 
R. In the following we shall neglect the distinction between S{t\R) and S{t) 
by assuming A <^r/g. 

Now, if the transformation a is such that (2.10) can be solved, at least in 
perturbation theory, denoting the solution by m(t) = h{t — to, mo) we have 

-1/ ,,dh~^{t — tQ,m) 



Pm(m, t) = Pmo(h (t - to, m)) 



dm 

(2.14) 

Po{a{h-\t - to, m))\R)^^^^^^t:^^a\h-\t - to, m)) 

dm 

which formally solves the problem since it gives an explicit expression for 
Pmijn-i t) in terms of the known quantities Po,a and h. 
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In order to see how this works in practice let us go back to the equation 
of motion (2.10) for m{t) and let us look for the transformation a which 
simplifies as much as possible the behaviour of m{t). The first attempt is for 
an outright linearization of (2.10). If this was not possible then, as stated 
above, a ought to be such that (2.10) can be solved in perturbation theory. 
However, in this case linearization can be achieved. Imposing 



a'{m) 

one finds 



[r — ga (m)] = rm (2-15) 



TTl 

4> = <m) = (2.16) 

[1 + {g/r)m^\ 2 

and 

m{t) = h{t — to,mo) = TTTiQ. (2-17) 

Indeed, we have that while m(i) grows exponentially (f){t) eventually satu- 
rates via (2.16) to the final equilibrium value 0eg = i-^r/gi. Namely, the 
transformation a accounts for the nonperturbative features of the problem. 

Putting together (2. 14), (2. 16) and (2.17) we have the exact solution of 
the problem in terms of the auxiliary variable m{t). The motivation for 
going to this form of the solution is that in more complicated cases where 
h and therefore Pm{m,t) cannot be explicitely obtained, the consideration 
that the auxiliary variable m(t) should not be much affected by the nonlinear 
nature of the problem authorizes to attempt a gaussian ansatz for PmiiTi, t). 
This will be the crucial step of the GAP approximation in the phase ordering 
problem. The difficulty with an ansatz however is that it may not be possible 
to control the corrections to it. In any case it should be clear that a gaussian 
ansatz does not amount to an overall linearization of the problem, since in 
(2.12) the ansatz amounts to use a gaussian form for P^{m,t) while the 
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nonlinearity remains through the exphcit factor o'^(m). To be more specific, 
since Pq is gaussian it is evident from (2.14) that Pm{m,t) is gaussian if a 
and h are hnear. Thus, should it be possible to find an expansion parameter 
A such that a and h become linear for A ^ 0, the gaussian approximation 
amounts to take this limit inside Pm{m,t) in (2.12) but not in the explicit 
factor a'^{m). 

Let us see how this works in the one particle context. Since in this case 
h{t — to, m) is already linear, in order to make Pmij^i t) gaussian we need to 
linearize only a in (2.14). From (2.16) we may write 

a\r~'m) = - '-j—f^^f^. (2.18) 

r [1 + [g/r)T ^m^\ 

and using this in (2.14) we obtain 

Pmim,t) = P^\m,t)K{T-'m,g/r) (2.19) 

with 
and 

expjg , .J 
K{r''m,g/r) = ^^"^^^''^^ '"''[ . (2.21) 

[1 + {g/r)r '^■m?\ 2 

Thus, in this case it is possible to identify the nonlinear coupling g with the 
expansion parameter A which generates gaussian statistics for m in the limit 
A ^ 0. Then, following the previous discussion, the lowest order is obtained 
by letting (7 — > in (2.19) but not in the explicit factor a'^{m). From (2.12) 
then we get 

dmP^^\m,t)a^{m) = -{1 - 0f^e^[l - erf{^)]} (2.22) 
-00 g 
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with y = r/(2Ar^5f) which gives S'^'^^t) ~ Ar^ at short time as in (2.7) and 
S^'^^t) ~ {r / g)[l — ^y7fy] for long time. The quahtative effect of the saturation 
is correctly reproduced, altough there is a quantitative discrepancy with (2.8) 
in the law of approach to equilibrium. In conclusion, in the one particle case 
the gaussian approximation can be identified with the zero order step in 
a systematic development where higher order corrections are generated by 
expanding K{T~^m, g/r) in powers of g. 



3 - Phase Ordering Dynamics 

Let us now turn to the field theory case. The phase ordering dynamics 
following the quench from high temperature to zero temperature of a system 
with a non conserved order parameter is described by the equation of motion 

^^^V'<f>{x,t)-V'{<P{x,t)) (3.1) 

with a gaussian initial state which generalizes (2.2) 

^'o[0o(^)] = ^e-^/<^''-*o(^) (3.2) 

and where V{(f)) is a potential of the double well type. 

Again, due to the symmetry of the problem the average order parameter 
vanishes identically < 0(x, t) >= and the observable of interest is the equal 
time correlation function 

G{u,t) =< (p{xi,t)4>{x2,t) >= j d(pid4>2P{4>i,Xit;(p2,X2t)4>i4>2 (3.3) 



or the structure factor 



C{k,t) = J (fxe^^-^G{u,t) (3.4) 
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where u = xi— X2- In (3.3) P(0i, Xit; 02, X2t) is the joint probabihty density 
that (p{x,t) takes the value 0i at the space-time point {xi,t) and the value 
02 at the space-time point {x2,t). 

It has been well established^'^, both from experiment and numerical sim- 
ulations, that in the late stage of the dynamics these quantities obey scaling 

G{u,t)r. f{u/L{t)) (3.5) 

C{k,t) ^ L\t)g{kL{t)) (3.6) 

where L{t) is the basic length in the problem which is related to the average 
size of domains and obeys the growth law L{t) ~ i^/^, while f{x) and g{x) 
are scahng functions. The origin of scahng is that in the late stage the order 
parameter reaches local equilibrium and forms domains of the ordered phases 
which evolve according to self-similar patterns. From the existence of sharp 
interfaces separating domains one can deduce^ the short distance behaviour 
of f{x) or the long wavelength behaviour of g{x) (Porod's law) 

f{u/L) = l-2u/L + ... foru/L«l (3.7) 

g{kL) ~ (A;L)-('^+^) for kL » 1 (3.8) 
as well as the saturation law"^ of the order parameter 

S{t) = G{u = 0, t) = <P%[1 - + 0{L'')] (3.9) 

where (peg is the value of the order parameter in the final equilibrium state. 
Eq.s from (3.5) to (3.9) contain the minimal phenomenological information 
that a theory of phase ordering dynamics should account for. 

At this point it is important to emphasize that the scaling behaviour 
described above applies to the late stage of the process where domains are 
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close to saturation and grow through the motion of the interfaces. This 
stage of the dynamics is dominated by the nonUnear nature of the problem 
and much as in the one particle case it cannot be obtained through any 
straightforward perturbation expansion. The great difference with the one 
particle case is that Eq. (3.1) cannot be solved for any realistic potential. 
Therefore, in order to make analytic progress, we turn to the generalization 
of the auxiliary variable method. 



4 - Auxiliary Field Method 

Following the idea illustrated above we now introduce an auxiliary field 
through a local nonlinear transformation 

(f){x,t) ^ a{m{x,t)) (4.1) 

which in general is defined through a relation involving the potential 

K[a{m)] - V'{a). (4.2) 

We note that such a transformation cannot be a linearising transformation 
as it was in the one particle case. In fact in that case m{x, t) ought to be the 
free field and the relation between the free field and the interacting field is 
certainly nonlocal, as it can be easily seen generating the formal solution of 
(3.1) by iteration. Thus, the transformation (4.1) is introduced in order to 
take care at least of the gross nonhnear effect which is the saturation of the 
order parameter to the finite final equilibrium value ^e?, leaving the rest, pos- 
sibly, to perturbation theory. Accordingly, for large time the transformation 
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must go over to the form 

a{m{x,t)) — (peg sign{m{x,t)). (4-3) 

The equation of motion of the auxihary field is obtained from (3.1) 

^ = V'm + l[a"(m)(Vm)' - V'{a)] (4.4) 

with the transformed initial condition 

Pmoi'moix)] = Po[(T(mo(f))]J(0o,"^o) (4.5) 

where J(0o, mo) is the Jacobian of the transformation (4.1) at the initial time. 
Representing the solution of (4.4) as a functional of the initial configuration 
labeled by x and t 

m{x, t) — h{x, t — to] [mo{x')]) (4-6) 

the probability of a configuration [^^(x)] at the time t can be obtained in 
terms of the initial probability density (3.2) 

Pm[m{x),t] = Pmo[h~^{x,t- to, [m(f )] )] J(mo, m) 

(4.7) 

Po[a{h~^{x,t - to, [m{x')]))]J{4>o,mo)J{mo,m) 

where is the inverse of (4.6) and J(mo, m) is the Jacobian of this transfor- 
mation. The above result is the analogue of (2.14) and specifies the statistics 
of the auxiliary field m{x, t) in terms of a, h and the statistical properties of 
the initial condition. 

Neglecting for simplicity considerations pertaining to the restriction of 
averages to domains of configurations where (4.1) is invertible, the correlation 
function (3.3) may be rewritten as 

G{u,t) = J dmidm2Pm{mi,xit;'m2,X2t)a{mi)a{m2) (4.8) 
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where the joint probabihty of m is related to (4.7) by 

J.i; m„ = / d[m(£)lP„[m(J), " m(£,))^(".. " 

(4.9) 

The above form (4.8) for the correlation function makes a progress over (3.3) 
if the joint probability of m is available. This requires that the transformation 
a is such that Eq.(4.4) for m is soluble. Short of this, as explained in section 
2, one resorts to the GAF approximation through the linearization of a and 
h inside Pq. 

Let us now review the predictions of the GAF approximations. If m(x, t) 
is gaussian the probability densities are of the form 

P^\mi,xit;m2,X2t) = ^exp{-— ^ [ml + ml - 2jmim2]} 

Zm 2(1 - T)So{t) 

(4.10) 

and 

P^\rn,xt) = ^ exp{--|^} (4.11) 

with 

'S'o(^) =< m^{x,t) >o , Go{u,t) =< m{xi,t)m{x2,t) >q 

(4.12) 



7 = 7(«,t) = ^|^ ' ^- = 2vr^o(t)Vl^ 

and where < • >o denotes averages with respect to ■ Hence, for the 
fluctuations of the order parameter one has 

Sit) = J dmPj^\m,x,t)a\m) (4.13) 

and for the scaling function 

/2 
dmidm2P^\m\,Xit\m2.,X2t)sign{mi)sign{m2) = —sin~^{^). 
TT 

(4.14) 
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Within this approach the problem is reduced to the computation of Gq{u, t). 
For h to be hnear Eq. (4.4) must be of the form 

dm o . X . ^ 

— = V^m + a(t)m (4.15) 
at 

where a{t) is some function of time to be determined. Upon hnearizing a the 
initial probability density (4.5) becomes 

Pmoimoi^)] - Po[cmo{x)] (4.16) 

where c is a constant. Solving (4.15) by Fourier transform and averaging 
over initial conditions with (4.16) one finds 

Co{k, t) = So{t)L''{t)go{kLit)) (4.17) 

Goiu,t) = Soit)-fi\u\/L{t)) (4.18) 

with 

' L{t) = 
go{kL) = exp(-2(A;L)2) 

So{t) = ^,exp{2b{t)) 
^ h{t)^!ldt'a{t'). 

Inserting the above expression for 7 in (4.14) one obtains the Ohta-Jasnow- 
Kawasaki^ result for the scaling function which correctly reproduces the be- 
haviours (3.7) (3.8). It is then matter of studying the behaviour of S{t) and 
for this we must go over to the specific implementations of the method. 

On site linearization 
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Making a direct extension to the field theory case of the procedure adopted 
for one particle, let us look for a transformation a which linearizes the on 
site potential in (4.4) 

- = rm (4.20) 

a'[m) 

where r is a constant. This yields 

— = V^m + rm- Q{m) {Vmf (4.21) 



where 



(7'(m) 

With the double well potential of the form 

^ x2 , 9 



= (4.22) 



V(<t>) = -^<P' + 1<I>' (4.23) 



(4.20) reduces to (2.15) yielding as in (2.16) 



m 



[1 + {g/r)Tm?]^ 



and 



Qim) = 3^^ , (4.25) 

However, contrary to what happens for one particle, even though the trans- 
formation (4.24) manages to account for the saturation of the order param- 
eter, yet it is not sufficient to linearize the equation of motion. This is 
done by introducing an approximation, which is optimized by the mean field 
prescription^^ 

g(m)(Vm)2 ^ 3^ < (Vm)2 >< ^ „ > (4.26) 

where averages must be computed self-consistenly. Note that although (4.26) 
yields the best linear approximation to the equation of motion, it remains an 
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uncontrolled approximation since no small parameter emerges which allows 
to compute corrections to it. According to the general discussion made above 
the implementation the GAF approximation requires, besides the lineariza- 
tion of the equation of motion, also the linearization of a. Setting a{m) — m 
in (4.24), the initial condition is given by (4.16) with c—\. 

With (4.26) the equation of motion is of the form (4.15) with 

a(i)=r-3^Z}o(t)|^ (4.27) 

where 

Do(i) (Vm)2 >= ik^Co{k,t). (4.28) 

Jk 

Next, using (4.13) 



and making the assumption to be verified a posteriori that So{t) grows with 
time, asymptotically we have 



S{t) = -{1 - J^r^ + Oi-^^}. (4.30) 
g V 2gSo{t) So{t) 



Inserting into (4.27), to dominant order we get 



K^)=r-3^ = r + 0(O (4.31) 

which gives b{t) = rt. Next, using (4.19) we find 

SM = A^25^ (4.32) 

which is consistent with the assumption made about So{t). Finally, inserting 
the above result into (4.30) we obtain that S{t) saturates exponentially fast 
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to the equilibrium value 4>lg = r/g, rather than according to a power law as 
expected from (3.9). 



KYG-theory 

The behaviour of S{t) obtained above is what one finds resumming the sin- 
gular perturbation series of Kawasaki, Yalabik and Gunton^^ (KYG). The 
KYG theory is contained in the above treatment as a particular case. If in 
addition to the mean field approximation one makes also an expansion in the 
nonlinear coupling g, to lowest order Q{m) = and the equation of motion 
becomes 



namely the auxiliary field coincides with the free field. The transformation 
(4.24) together with (4.33) corresponds exactly to the KYG theory, which 
therefore in the present context amounts to the statement that all the im- 
portant nonlinear features of the problem are adequately taken care of by 
the transformation (4.24). 

BH-theory 

If in (4.26) we keep the first order in g the equation of motion becomes 



which is of the type of the equation obtained by Bray and Humayun^ (BH) 
starting from an ad hoc potential and which leads to the correct behaviour 
for S{t). In fact in this case (4.27) reduces to 




(4.33) 



dm 
'dt 



V^m + [r-3- < {Vmf >]m 



(4.34) 




(4.35) 
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and setting to zero the left hand side for large time 



(4.36) 



Inserting this result in (4.30) the behaviour (3.9) of S{t) is recovered. This 
is due to the cancellation of So{t) in the denominator of (4.27), which occurs 
only in first order in g. Notice that from the above result for So{t) and the 
definition (4.19) one obtains a{t) ~ {d + 2)/4t which coincides with the form 
for a{t) introduced by Oono and Puri^^ in their improvement of the OJK 
theory. 

In summary, the GAF approximation obtained via the linearization of the 
on site potential i) does not describe correctly the saturation law of the order 
parameter, except for the very special case where the BH-theory applies, 
and ii) it is an uncontrolled approximation since there is not a systematic 
expansion scheme within which it plays the role of the zero order theory. 

Mazenko transformation 

Let us now go to a different way of introducing the auxiliary field due to 
Mazenko^ where equation (4.2) is chosen in such a way that a{m) reproduces 
the profile of the static interface 



In this case m{x, t) has the physical interpretation of the distance to the 
nearest interface. Using (4.37) in (4.4) we obtain 



a"{m) = V'{a). 



(4.37) 



dm 
'dt 



V^m+ (1 - {Vmf)Q{m) 



(4.38) 
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where Q{m) is still given by (4.22). Note that since V'{a) is an odd function 
from (4.37) and (4.22) follows that also Q{m) is an odd function. Thus, the 
mean field linearization of (4.38) yields 



where H{t) is some function of time whose explicit form is not important. 
Hence Eq. (4.27) now gives 



which, apart for the overall factor H{t), is identical to (4.35) and therefore 
leads to the same result (4.36) for So{t) which yields the correct behaviour of 
S{t). Comparing (4.34) with (4.39) we see that the BH-theory is a particular 
case arising with H{t) constant. Thus, the GAF approximation obtained 
within the static interface approach yields correct results, but for the same 
reasons pointed out above it remains an uncontrolled approximation. 

Generalization of the transformation 

We end up this section by considering a generalization of the transformation 
obtained by allowing for an explicit time dependence. The idea is to see if in 
so doing one may get closer to the linearization of the equation for m{x,t), 
as suggested by recent work of Puri and Bray^''. Replacing (4.1) by 



dm 
'dt 



V^m + [1- < (Vm)^ >]H{t)m 



(4.39) 




(4.40) 



0(f , t) 



a{t, m{x, t)) 



(4.41) 



the equation for m becomes 




da /dm dm^ 



[7^(Vm)^ 



da^ 
'dt 



V'{a)]. 



(4.42) 
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Let us then determine the exphcit dependence of cr on t by imposing 

% = -V'{cj) (4.43) 

which yields 

= V^m - Q{t, m)(Vm)2 (4.44) 

with 

da I am 

Eq. (4.43) is nothing but the one particle equation of motion which, using 
the potential (4.23), yields the solution (2.4) i.e. 

Tg({). Ill) 

a(t,m) = 4.46 

[l + %2(o,^)(^2_ip 

and 

Ir r J 

where the sigma's on the right hand side stand for a{Q, m) and the primes 
denote derivatives with respect to m. 

Imposing cr'V^ — ?>a{a'Y = we find (7(0, m) = ±(m)~^/^ and inserting 
into (4.47) and (4.46) eventually we have 

2 

T TTX 1 

1 + {g/r)m{T^ — 1) 

and 

dt 2 [r-2m- ^(t-2 - 1)]^ ^ ^ ^ 

Indeed, the equation of motion for m is "almost" linear since the nonlinear 
term vanishes exponentially fast, but the scheme it is not of much use in 
generating a GAP approximation, since (4.48) cannot be linearized. 
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5 - Vector Fields 



Let us now consider the case of a vector order parameter with A'"-components 
(j){x) = (01 (x), 4>n{x)). In this case a systematic expansion scheme about 
the GAF approximation can be outhned, although its practical implementa- 
tion remains to be explored. 

Phenomenological expectations in this case are a power law tail in the 
scaling function of the structure factor^'^ 

gix) ~ x-(<^+^) (5.1) 

which generalizes Porod's law and the saturation law^ 

S{t) = <PUl-j^^+0{L--')] (5.2) 

in place of (3.9). Considering the equation of motion 

= VVa(f , t) - ^^(fe t)) (5.3) 

with the potential 

Vi^) = + i^i^r (5.4) 
the auxiliary field m{x, t) is introduced by generalizing to the vector case the 
transformation (4.24) 

(^a(rn) = — r (5-5) 

which yields the equation of motion for m 

■ma + rnia - vm^j 



(5.6) 



[2Vma • J2'y{'m-yVm^) - ^ma(m-^Vm^)^] 
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For the equal time correlation function G{u, t) =< 0q(5i, t)0a(a;2, t) >, which 
is independent of a due to the rotational symmetry of the potential, we have 

G{u,t) = y" (imi(im2-Pm(^i) m2, a;2t)cra(mi)aa(m2) (5.7) 

where Pm, which is related to the initial probability density Pq through the 
analogues of (4.7) and (4.9), depends explicitely on A'^ through a and h. As 
previously stated Pm becomes gaussian upon linearizing a and h. We now 
show that this is achieved by taking the large- limit. The major difference 
with the scalar case than is that now there emerges A = as the natural 
parameter which yields the gaussian approximation in the limit A ^ 0. 

Taking the limit N ^ oo terms of the type J2a Qa in (5.5) and (5.6) are 
replaced by the average < Qa > yielding the linear equations 

CTairn) = — r (5-8) 

and 

^ = V'm« + [r - ^ < (Vm^)" >]m« (5.9) 

since < m^Vm^ > vanishes. Hence, as anticipated, in the large- limit 
the auxiliary field m is gaussian. Furthermore Eq. (5.9) is of the BH-type 
yielding (4.36) for So{t) =< ml{t) >. 

It is important to realize that the large- iV limit we are considering here 
is quite different from the usual large- A?^ limit^^ performed on the equation 
of motion (5.3) for 0. The latter one is recovered in the present context 
by taking the large- A?" limit, namely using (5.8), also in the explicit a^s 
appearing in (5.7) and eventually obtaining 

G(^?.^) = ^[l-^]exp(-^). (5.10) 
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Instead, according to the general structure of the GAF approximation which 
we have repeatedly illustrated above, A^ must be kept fixed to whatever value 
it has been originally specified in the explicit cTq 's in (5.7), while the N ^ oo 
limit is taken inside Pm. In so doing from (5.5) and (5.7) we obtain the Bray, 
Puri and Toyoki ^ (BPT) result for the scahng function 



u 

fij^) =< mixi, t) ■ m{x2, t) > 



(5.11) 



A^7 



whereB{x,y) is the beta function, F{a,b;c;z) the hypergeometric function 
and ^{u/L) is given by (4.19). Prom the above result follows the power law 
tail (5.1). Furthermore, from S{t) =< > we obtain 

, 1 f dm _^ , , 

and carrying out the integral 

Nr Nr m/o jvr N Nr 

where T{x, y) is the incomplete gamma function. Expanding up to first order 
in l/5'o we obtain 

f 4>Ul - T^^^\ for AT > 2 
S{t) - <^ (5.14) 
i '^ejl + ^] for = 2 

which yields the power law behaviour (5.2), contrary to the exponential sat- 
uration which one obtains in the BTP approach. 

These results show that the expected phenomenological behaviour is ob- 
tained at zero order within the 1/A^-expansion of the probability density of 
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the auxiliary field. In principle, systematic corrections could be obtained 
via the higher order terms in the l/A'^-expansion of Pm, although we do not 
expect that such a scheme of computation might be easily implemented in 
practice. It is worth pointing out that the scheme for the systematic im- 
provement of the GAF approximation for vector fields presented here is con- 
ceptually different from that proposed by BH in two respects: i) while we use 
the standard quartic potential (5.4) BH need to invoke an ad hoc potential 
which cannot even be written in closed form and ii) the expansion is made 
in where here is the number of components of the order parameter 
rather than the number of components of an additional internal color index. 
Finally, the comparison between (5.10) on one side and (5.11) (5.14) on the 
other clearly shows the difference between the standard l/A"-expansion and 
the one we have presented here. The most important point is that while 
there are no localized defects in lowest order in the usual l/A'-expansion 
since the correlation function (5.10) decays exponentially, the power law tail 
(5.1) implied by (5.11) shows that our reformulation of the l/A"-expansion 
describes defects in lowest order. 



6 - Concluding remarks 



In conclusion, in this paper we have analysed the sequence of steps which 
must be taken within the framework of a first principles theory in order 
to generate GAF approximations. The analysis has been restricted to sys- 
tems with non conserved order parameter. The idea was to look for the 
systematic expansion scheme which allows to control the corrections to the 
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GAF approximations. A project of this type is suggested by the physical 
motivation behind the introduction of the auxihary field. This being more 
smooth and less non linear than the order parameter field, hopefully should 
be tractable in perturbation theory. Our results are negative for the scalar 
case, in the sense that we are unable to come up with the expansion scheme 
within which the GAF approximation can be identified with the zero order 
approximation. It should be mentioned that there are indications^^'^ that the 
GAF approximation becomes exact in the limit of infinite space dimensional- 
ity, suggesting the 1 / rf-expansion as a possible systematic expansion scheme. 
This is an interesting line of research worth to be further investigated. 

The outlook is somewhat better in the case of a vector order parameter. 
In this case one can set up the theory in such a way that the large- A?^ limit 
yields the GAF approximation. Consequently one can expect that there 
exists a l/A"-expansion where corrections to the GAF approximation are 
generated systematically. Finally, approximations which go beyond the GAF 
approximation have been introduced recently by Mazenko^^. In future work 
we plan to look for the connection between that work and the point of view 
developed here. 
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